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ABSTRACT 



We present the results of a CO line survey of 30 galaxies at moderate redshift (z ~ 0.2-0.6), with the IRAM 30m telescope, with 
the goal to follow galaxy evolution and in particular the star formation efficiency (SFE) as defined by the ratio between far-infrared 
luminosity and molecular gas mass (Lfir/M(H2))- The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with 
Lpi R larger than 2.8 10 12 L G , experiencing starbursts; adopting a low ULIRG CO-to-H 2 conversion factor, their gas consumption time- 
scale is lower than 10 s yr. To date only very few CO observations exist in this redshift range that spans nearly 25% of the universe's 
age. Considerable evolution of the star formation rate is already observed during this period. 18 galaxies out of our sample of 30 are 
detected (of which 16 are new detections), corresponding to a detection rate of 60%. The average CO luminosity for the 18 galaxies 
detected is L' C o = 2 10'° K kms~' pc 2 , corresponding to an average H 2 mass of 1.6 10 10 M Q . The FIR luminosity correlates well 
with the CO luminosity, in agreement with the correlation found for low and high redshift ULIRGs. Although the conversion factor 
between CO luminosity and H 2 mass is uncertain, we find that the maximum amount of gas available for a single galaxy is quickly 
increasing as a function of redshift. Using the same conversion factor, the SFEs for z~0. 2-0.6 ULIRGs are found to be significantly 
higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected 
cosmic H 2 abundance and the cosmic star formation history. 

Key words. Galaxies: high redshift — Galaxies: ISM — Galaxies: starburst — Radio lines: Galaxies 



1. Introduction 

Ultra-Luminous Infra-Red Galaxies (ULIRGs) emit most of 
their energy in the far-infrared, and have far-infrared luminosi- 
ties L F ir > 10 12 L , (e.g. Sanders & Mirabel 1996, Veilleux etal. 
2009). Since they can be seen so far away, they allow us to ex- 
plore the evolution of star formation in the universe, and of the 
star formation efficiency (SFE) in particular, defined as the FIR 
luminosity to H2 mass ratio (e.g. Kennicutt 1998). Since the dis- 
covery of the first high redshift object in CO line emission (IRAS 
F10214+4724 at z=2.3, Brown & Vanden Bout 1991, Solomon 
etal. 1992), there has been a wealth of CO-line discoveries, a 
hundred objects are now detected at z>l, either from ULIRGs, 
or from LIRGs (Lfir > 10 11 L Q ). Some are amplified by grav- 
itational lensing (see the review by Solomon & Vanden Bout 
2005). They allow us to observe the interstellar medium of the 
galaxies, the CO excitation (e.g. Weiss etal. 2007) and estimate 
the amount of molecular gas present. Stars form from molecular 
gas, so it is important to infer the H2 mass in order to determine 
the SFE. At high redshift, many of these objects are quasars or 
radio-galaxies (due to their selection, e.g. Omont etal. 2003), 
however, their FIR emission is powered predominantly by star 
formation (e.g. Riechers etal. 2009, Wang etal. 2010). 
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Locally, our knowledge of the ULIRG phenomenon is 
more profound due to higher spatial resolution and sensitivity. 
Because the peak of the dust emission is progressively shifted 
from the FIR to the submm domain, the dust emission can be 
detected to high redshifts (negative K-correction, e.g. Blain & 
Longair 1996). The CO-line emission is less favoured, and CO 
lines are difficult to detect at high z although observing the high- 
J CO lines helps significantly in highly excited objects (Combes 
etal. 1999). To date, more than a hundred objects have been 
studied in detail locally. In the case of the ULIRGs, it was found 
that they are characterized by compact, nuclear starbursts (e.g. 
Downes & Solomon 1998), and it has been argued that a spe- 
cial CO-to-H2 conversion factor should be used, that is 5.75 
lower than the standard factor commonly used for Milky Way- 
like galaxies (Downes et al. 1993). In the present paper, we will 
adopt for ULIRGs the ratio a=0.8 (Solomon etal. 1997) be- 
tween M(H2) and L'co, expressed in units of M (K kms -1 
pc 2 ) 1 , and not the standard a-4.6. 

At intermediate redshifts, between 0.2 < z < 1, there is 
a dearth of CO-line detections. This is partly due to observa- 
tional difficulties. The most commonly used millimetric win- 
dow is the 3mm one, which is least affected by atmospheric 
opacity. Between 81 and 115 GHz, all redshifts can be observed 
with at least one line of the CO rotational ladder, except be- 
tween z=0.4 and 1 . The latter can be observed at 2mm (targeting 
the CO(2-l) line), but in less favorable atmospheric conditions. 
While between z=0.2 and z=0.4 the 3mm window can be used 
in the CO( 1 -0) transition, the K-correction is strongly reducing 
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its observable intensity (by a factor growing faster than (1+z) 4 , 
Combes etal. 1999). The redshift range 0.2 < z < 1 is impor- 
tant though, as it covers almost half of the age of the universe, 
and also the most dramatic change in star formation activity (e.g. 
Madau etal. 1998, Hopkins & Beacom 2006). In the universal 
star formation history, the most striking feature is the impressive 
drop between z=l and z=0 by at least an order of magnitude 
(Blainetal. 1999). 

Up to now, very little was known about the molecular gas 
content of galaxies at moderate redshift. The ULIRG sample 
of Solomon etal. (1997) contains 37 objects, but only 2 have 
Z > 0.2. Negative results were obtained in previous studies, 
conducted about lOyrs ago (Lo etal. 1999, Wilson & Combes 
1998), but the performances of the mm-instruments have dra- 
matically improved since then. Two more objects were detected 
by Geach et al. (2009), although more upper limits were also re- 
ported (Melchior & Combes 2008). To study star forming galax- 
ies in this period, and in particular to derive their star forma- 
tion efficiency, we have undertaken a CO-line search in the range 
0.2 < z < 0.6, almost unknown territory as far as molecular lines 
are concerned. We have selected a sample of 30 IR-luminous 
galaxies in this redshift range to check whether the variation 
of star forming activity is due to a variation in molecular gas 
content or star formation efficiency, or both. One of the objects 
(IRAS 11582+3020, hereafter G4) has already been mapped 
with the IRAM Plateau de Bure Interferometer (Combes et al. 
2006, paper I). The CO map showed spatially resolved emission 
on 30kpc scales and revealed a velocity gradient. It was con- 
cluded in that paper that not all the molecular gas is confined 
in a nuclear starburst, but that ~50% of it is extended on galac- 
tic scales (25-30kpc). In the present paper, we describe the CO 
survey carried out with the IRAM 30m telescope. The sample is 
described in Sect. [2] and the observations in Sect. [3] Results are 
presented in Sect. [4] and discussed in Sect. [5] 
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Fig. 2. The CO spectra of the detected galaxies. The zero veloc- 
ity scale corresponds to the optically determined redshift, listed 
in Table Q] Sources detected in CO(3-2) at 1mm wavelength are 
also shown. Some sources were detected in CO(1-0) but not in 
CO(3-2), as indicated in Table [2] The vertical scale is T m t in 
mK. The spectrum of G4 is already presented in Paper I. 
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Fig. 1. Definition of our sample. Among the 209 northern galax- 
ies (filled symbols) found in NED between 0.2 < z < 0.6 and de- 
tected at 60 //m by IRAS, we selected the most luminous ones 
(log Lfir/Lq > 12.45, as indicated by the horizontal line). By 
comparison, the ULIRGs in the sample of Solomon et al. (1997) 
are plotted as open triangles. The circles indicate detections, 
non-detections are marked by a cross. Sources that have neither 
a circle nor a cross could not be observed due to weather condi- 
tions. 



2. The sample 

The present-day sensitivity in the CO line restricted the sample 
to the brightest objects in the far infrared. We have selected all 
objects between 0.2 < z < 0.6 and DEC(2000) > -12° that are 
identified as galaxies in the literature, have spectroscopic red- 
shifts, and are detected at 60/mi (IRAS, ISO). This resulted in a 
total of 209 galaxies. Most of the galaxies (and in particular the 
brightest ones) have detailed photometry in the NIR bands, from 
the samples by Clements etal. (1996), Kim & Sanders (1998), 
Kim etal. (2002) and Stanford etal. (2000). The available sub- 
arcsec K-band images (from either IRTF or Keck telescopes, in 
the above references), reveal that about two-thirds of the objects 
are interacting galaxies. 

Out of the 209 galaxy sample, we selected the brightest ones, 
with log Lfir > 12.45. This leads to a sample of 36 objects, to 
be observed with the 30m telescope. We did not reobserve one 
galaxy detected by Solomon et al. (1997), nor 3C48, detected by 
Wink etal. (1997), although we include these 2 sources in our 
analysis. Due to weather conditions, only 28 sources were ob- 
served in the project. All their identification and coordinates are 
displayed in TableQ] Out of the 30 objects in the sample (includ- 
ing the two literature sources), 18 were detected, corresponding 
to a detection rate of 60%. Figure Q] displays the distribution of 
FIR luminosities with redshift. 

The far-infrared fluxes Ffir are computed as 1.26 x 10" 14 
(2.58 Seo+Sioo) W nr 2 (Sanders & Mirabel 1996). The far- 
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Fig. 3. Same as Fig.[2]for the following galaxies. 

infrared luminosity is then Lfir = 4 ji D 2 L CC Ffir, where 
D/, is the luminosity distance, and CC the color correction, 
CC=1.42 (e.g. Sanders & Mirabel 1996). The FIR-to-radio ra- 
tio <7=log([F FIR /(3.75 10 12 Hz)]/[f v (1.4 GHz)]) has been com- 
puted for sources where radio data were available; the radio 
fluxes are listed in Table [2] Excluding the radio galaxies 3C48 
and 3C345, the average is q-2.3, typical for ULIRGs (Sanders 
& Mirabel 1996). The star formation rates of all sample galaxies 
are above 480 M Q yr -1 , estimated from the infrared luminosity 
(e.g. Kennicutt 1998). 

In this article, we adopt a standard flat cosmological model, 
with A = 0.73, and a Hubble constant of 71kms~ 1 Mpc~ 1 
(Hinshawetal. 2009). 

3. Observations 

The observations were carried out with the IRAM 30m telescope 
at Pico Veleta, Spain, between January 2005 and January 2006. 
Most of the galaxies, with redshifts between 0.2 and 0.39, could 
be observed simultaneously at 3mm in CO(1-0) and at 1mm in 
CO(3-2), except the two lowest redshift sources (G4 and G5) 
where only observations in the 3mm band were possible. In any 
case, the weather prevented us sometimes from taking useful 
data at 1mm. For the highest redshift sources (G19, G20, G23, 
G27 and G28), only the CO(2-l) line was observed in the 2mm 
band. 

The SIS receivers were tuned in single sideband mode to the 
redshifted frequencies of the various CO lines. The observations 
were carried out in wobbler switching mode, with reference po- 
sitions offset by 3' in azimuth. We used the 1 MHz back-ends 
with an effective total bandwidth of 512 MHz at 3 mm (provid- 
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Fig. 4. Same as Fig.|2]for the remaining galaxies. 

ing a ~ 1500 kms" 1 range) and the 4 MHz filterbanks with an 
effective total bandwidth of 1024 MHz at 1 mm. 

We spent 2-A hours on each galaxy, resulting in a relatively 
homogeneous noise level of 1-2 mK per 30 km s _I channel for 
all sources. The system temperatures ranged between 120 and 
250 K at 3 mm, between 220 and 300 K at 2 mm, and between 
300 and 500 K at 1.2 mm, in T^. The pointing was regularly 
checked on continuum sources and yielded an accuracy of 3" 
rms. The temperature scale used is in main beam temperature 
r m b- At 3mm, 2mm and 1mm, the telescope half-power beam 
width is 27", 17" and 10" respectively. The main-beam efficien- 
cies are rj mb = T^/T m b=0.S5, 0.70 and 0.64, respectively, and 
S /T mb =4.8 Jy/K for all bands. 

Each spectrum was summed and reduced using linear base- 
lines, and then binned to 50 - 60 km s~' channels for the plots. 

4. Results 

4.1. CO detection in z=0.2-0.6 ULIRGs 

All spectra for CO detections are displayed in Figures [2] [3] and [4] 
(except G4 reported in paper I). The non-detections are reported 
in Table [2] Integrated upper limits are computed at 3cr, assum- 
ing a common line-width of 300 kms -1 , and getting the rms of 
the signal over 300 km s~ 1 . Lines are assumed detected when the 
integrated signal is larger than 3cr. Gaussian fits then yielded the 
central velocities, velocity FWHMs and integrated fluxes listed 
in Table [2] 

As already noticed in Sect. |2l very few objects were pre- 
viously detected in CO in this redshift range. We include in 
our analysis, and in Table [2] two additional ULIRGs (G29 and 
G30) that satisfy our sample criteria. In the discussion, we also 
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Table 2. Observed line parameters 



Galaxy 


Line 


y obs 


S(CO)° 


V b 


AVpWHM 


L'co/10 10 


Seo 


100 


log LpiR 


F(1.4GHz) rf 






[GHz] 


[Jy kms - '] 


[kms -1 ] 


[kms 




[K kms - ' pc 2 ] 




rjvi 

L J j J 




taJvl 

L±±±j j j 


Gl 


CO(l-O) 


88.534 


7.5 ± 1.4 


-29. ± 33. 


398. ± 


98. 


3.58 


0.68 


0.77 


12.54 




Gl 


CO(3-2) 


265.588 


6.2 ± 2.0 


-38. ±21. 


148. ± 62. 


0.33 


0.68 


0.77 


12.54 




G2 


CO(l-O) 


83.228 


< 1.2 








< 0.9 


0.67 


1.13 


12.84 


4.09 


G2 


CO(3-2) 


249.672 


<3.2 








<0.3 


0.67 


1.13 


12.84 


4.09 


G3 


CO(l-O) 


92.513 


2.3 ± 0.5 


-168. ±41. 


314. ± 


66. 


0.73 


1.18 


1.55 


12.59 


2.16 


G3 


CO(3-2) 


277.525 


9.6 ± 0.9 


-185. ± 37. 


498. ± 


47. 


0.33 


1.18 


1.55 


12.59 


2.16 


G4 


CO(l-O) 


94.253 


6.5 ± 0.7 


248. ± 24. 


434. ± 


52. 


1.67 


1.23 


1.52 


12.50 


3.09 


G5 


CO(l-O) 


94.718 


3.4 + 0.6 


178. ± 20. 


240. ± 


44. 


0.83 


1.36 


1.54 


12.51 


25.5 


G6 


CO(l-O) 


90.409 


< 1.7 








<0.7 


0.72 


0.69 


12.44 


3.0 


G7 


CO(l-O) 


84.758 


< 1.7 








< 1.2 


1.18 


1.20 


12.94 




G7 


CO(3-2) 


254.262 


<4.5 








<0.3 


1.18 


1.20 


12.94 




G8 


CO(l-O) 


93.186 


3.4 ± 0.7 


-135. ±32. 


301. ± 


82. 


1.0 


0.99 


1.43 


12.49 


5.96 


G9 


CO(l-O) 


92.439 


2.1 +0.5 


469. ± 30. 


290. ± 


72. 


0.67 


0.94 


1.10 


12.48 


3.85 


G10 


CO(l-O) 


91.850 


< 1.7 








<0.6 


0.82 


1.04 


12.46 


7.44 


Gil 


CO(l-O) 


89.985 


4.5 ± 0.6 


85. ± 30. 


428. ± 


66. 


1.86 


0.84 


1.15 


12.58 


2.77 


G12 


CO(l-O) 


92.811 


7.7 + 0.8 


-5. ± 21. 


430. ± 


44. 


2.33 


1.48 


1.99 


12.68 


5.67 


G13 


CO(l-O) 


88.875 


< 1.7 








<0.8 


1.24 


2.13 


12.84 


95. 


G14 


CO(l-O) 


92.811 


2.4 ± 0.7 


-10. ± 70. 


420. ± 


140. 


0.73 


1.06 


1.43 


12.53 




G15 


CO(l-O) 


88.330 


6.9 ± 1.1 


200. ± 34. 


433. ± 


73. 


3.38 


1.22 


1.27 


12.79 




G16 


CO(l-O) 


86.932 


< 1.2 








<0.7 


0.58 


1.14 


12.63 


40.3 


G17 


CO(l-O) 


82.337 


8.3 ± 0.9 


478. ± 19. 


363. ± 


46. 


7.14 


0.32 


0.64 


12.59 




G17 


CO(3-2) 


246.997 


12.9 ± 2.0 


500. ± 22. 


303. ± 


56. 


1.24 


0.32 


0.64 


12.59 




G18 


CO(l-O) 


85.008 


7.6 ± 0.7 


253. ± 32. 


651. ± 


64. 


5.18 


0.61 


1.22 


12.75 


2.3 


G18 


CO(3-2) 


255.012 


<4.5 








<0.3 


0.61 


1.22 


12.75 


2.3 


G19 


CO(2-l) 


159.874 


1.5+0.4 


-28. ± 30. 


225. ± 


63. 


0.40 


0.53 


<0.44 


12.78 


6.87 


G20 


CO(2-l) 


154.723 


< 1.4 








<0.5 


0.23 


0.57 


12.70 


8.82 


G21 


CO(l-O) 


88.670 


< 1.7 








<0.8 


0.71 


0.76 


12.54 


7.78 


G21 


CO(3-2) 


265.997 


<4.5 








< 0.2 


0.71 


0.76 


12.54 


7.78 


G22 


CO(l-O) 


87.129 


1.5+0.3 


-41. ± 26. 


196. + 


44. 


0.82 


0.62 


0.99 


12.62 


1.59 


G22 


CO(3-2) 


261.373 


<2.3 








<0.1 


0.62 


0.99 


12.62 


1.59 


G23 


CO(2-l) 


159.763 


4.7 ± 0.7 


91. ±25. 


370. ± 


60. 


1.27 


0.23 


0.62 


12.61 


1.76 


G24 


CO(l-O) 


85.134 


<0.8 








< 0.6 


0.52 


1.83 


12.80 




G24 


CO(3-2) 


255.388 


<4.5 








<0.3 


0.52 


1.83 


12.80 






V. V V l 1 — U 1 


/ .OJ7 










<- O f. 




i nf, 

1 ,UU 






G25 


CO(3-2) 


263.564 


<9.0 








<0.5 


0.40 


1.06 


12.49 




G26 


CO(l-O) 


86.410 


< 1.2 








<0.7 


0.46 


0.83 


12.54 


1.36 


G26 


CO(3-2) 


259.217 


<4.5 








<0.3 


0.46 


0.83 


12.54 


1.36 


G27 


CO(2-l) 


146.373 


2.1 ±0.6 


-126. ± 34. 


231. ± 


90. 


0.95 


0.43 


<0.94 


13.11 


2.7 


G28 


CO(2-l) 


144.719 


< 1.9 








<0.9 


0.60 


1.26 


13.28 


7000 


029' 


CO(l-O) 


84.170 


1.9 + 0.3 


-6. ± 10. 


270. ± 


20. 


1.42 


0.74 


0.83 


12.78 


16000 


G3f/ 


CO(l-O) 


91.123 


3.8 + 0.4 


185. ± 20. 


270. ± 


30. 


1.40 


1.45 


1.82 


12.75 


5.49 



Quoted errors are statistical errors from Gaussian fits. The systematic calibration uncertainty is 10%. 
" The upper limits are at 3tr with an assumed AV = 300 km s -1 . * The velocity is relative to the optical redshift given in TableQ] 
The 60 and lOOyurn fluxes are from NED (http://nedwww.ipac.caltech.edu/) 
d From the FIRST catalog (http://sundog.stsci.edu/). Errors are typically 0.14 mjy 
e From Wink etal. (1997). 3 From Solomon etal. (1997). 



added the two galaxies on the outskirts of the cluster CI 0024+16 
(Geach et al. 2009); they are not ULIRGs, but it is interesting to 
compare star formation efficiencies for all CO-detected objects 
in this redshift range. 

The detection rate of 60% in our sample, down to a sen- 
sitivity limit of ~ 1.5 Jykms~', must be considered a lower 
limit. Indeed, the available velocity range of the receivers (about 
1500 km s -1 ) could have missed some sources if the optical red- 
shift was not accurate enough. Some galaxies show a signifi- 
cant velocity offset (e.g. G17) as can be seen in Table [2] and 
the figures. Some of the profiles may have a double-horn shape 
as G18, but most do not, given our spectral resolution and sen- 
sitivity. The line-widths detected are compatible with massive 
galaxies at random inclinations. Their average is AVfwhm= 348 



km s _1 , very similar to the value for local ULIRGs of 302 km s _1 
(Solomon etal. 1997). In comparison, the submillimeter galax- 
ies have much broader widths, 655 kms -1 on average (Greve 
2005). Given the angular distance of the sources (average value 
1000 Mpc), our beam subtends between 50 and lOOkpc, and all 
galaxies can be considered unresolved, at least as far as their 
molecular component is concerned. 



4.2. CO luminosity and H 2 mass 

To derive the total H2 mass, we first compute the CO luminosity 
through integrating the CO intensity over the velocity profile. 
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Table 3. Molecular gas mass and star formation efficiency 



u 


Source 


KA(2U0U ) 


UnL(zUOU) 


z 


O 


M(ri2) 


orb 


i,i 


M d 


Type 














10' M Q 


L /M 


K 


10 s M Q 




Gl 


IRAS 00302+3625 


00:32:57.6 


+36:41:56 


0.3023 






G2 


IRAS 08081+2611 


08:11:14.4 


+26:02:17 


0.3850 


Gl 


28.6 


120. 


51.6 


0.8 




G3 


IRAS 10091+4704 


10:12:16.7 


+46:49:43 


0.2460 


G2 


<7.5 


>910. 


47.0 


2.6 




G4 


IRAS 11582+3020 


12:00:46.8 


+30:04:15 


0.2230 


G3 


5.8 


667. 


46.5 


1.6 


L,Int 


G5 


"J12054771+1651085 


12:05:47.7 


+ 16:51:08 


0.2170 


G4 


13.4 


238. 


46.7 


1.2 


L 


G6 


6 J 1307006+233805 


13:07:00.6 


+23:38:05 


0.2750 


G5 


6.7 


483. 


48.1 


1.1 


L, Int 


G7 


"J13301520+3346293 


13:30:15.2 


+33:46:29 


0.3600 


G6 


<5.3 


>521. 


54.4 


0.5 


Sl.Q 


G8 


IRAS 13352+6402 


13:36:50.7 


+63:47:03 


0.2366 


G7 


<9.4 


>928. 


56.5 


1.2 




G9 


IRAS 13379+3339 


13:40:14.4 


+33:24:45 


0.2473 


G8 


7.9 


391. 


44.5 


1.6 


Pair 


G10 


IRAS 13447+2833 


13:47:05.5 


+28:18:05 


0.2551 


G9 


5.3 


565. 


48.7 


0.9 


S2, Int 


Gil 


IRAS 15298+6319 


15:30:41.1 


+63:09:40 


0.2810 


G10 


<4.6 


>637. 


47.4 


1.0 




G12 


IRAS 16300+1558 


16:32:21.4 


+ 15:51:45 


0.2417 


Gil 


14.9 


255. 


46.9 


1.4 


Int 


G13 


[HB89] 1821+643 


18:21:57.3 


+64:20:36 


0.2970 


G12 


18.6 


255. 


46.0 


2.0 


L, Int 


G14 


IRAS 20551+2441 


20:57:19.7 


+24:53:37 


0.2425 


G13 


<6.3 


>1110. 


43.8 


4.0 


S1,Q 


G15 


IRAS 23113+0314 


23:13:54.3 


+03:30:58 


0.3053 


G14 


5.8 


584. 


45.7 


1.5 




G16 


IRAS 01506+2554 


01:53:28.3 


+26:09:40 


0.3264 


G15 


27.0 


227. 


53.5 


1.1 




G17 


IRAS F02 11 5+0226 


02:14:10.3 


+02:40:00 


0.4000 


G16 


<5.3 


>802. 


42.7 


2.9 




G18 


IRAS 07449+3350 


07:48:10.6 


+33:43:27 


0.3560 


G17 


57.1 


67. 


44.9 


2.0 




G19 


6 J09 13454+405628 


09:13:45.4 


+40:56:28 


0.4420 


G18 


41.4 


134. 


43.5 


3.4 




G20 


IRAS F10156+3705 


10:18:34.5 


+36:49:52 


0.4900 


G19 


3.2 


1875. 


>60.0 


<0.5 


Q2, Int 


G21 


IRAS 12514+1027 


12:54:00.8 


+ 10:11:12 


0.3000 


G20 


<3.8 


>1305. 


44.6 


2.8 




G22 


[HB89] 1402+436 


14:04:38.8 


+43:27:07 


0.3233 


G21 


<6.4 


>542. 


52.9 


0.7 


S2, Int 


G23 


C J145658.42+333710.1 


14:56:58.4 


+33:37:10 


0.4430 


G22 


6.6 


626. 


45.9 


1.8 


S1,Q, Int 


G24 


IRAS 19104+8436 


19:01:44.5 


+84:41:25 


0.3544 


G23 


10.2 


400. 


42.1 


3.3 




G25 


IRAS F00415-0737 


00:44:05.6 


-07:21:13 


0.3120 


G24 


<4.5 


>1379. 


36.6 


11.9 


S1,Q 


G26 


'7020412.43-00535 1.4 


02:04:12.4 


-00:53:51 


0.3343 


G25 


<4.9 


>631. 


38.5 


4.0 




G27 


IRAS F00235+ 1024 


00:26:06.5 


+ 10:41:32 


0.5750 


G26 


<5.6 


>615. 


44.2 


1.9 




G28 


3C345 


16:42:58.8 


+39:48:37 


0.5928 


G27 


7.6 


1699. 


>49.0 


<4.2 




G29 


3C48 


01:37:41.3 


+33:09:35 


0.3695 


G28 


<7.3 


>2615. 


50.3 


5.3 


Q 


G30 


f J140931.25+051 131.2 


14:09:31.2 


+05:11:31 


0.2644 


G29 


11.4 


523. 


54.5 


1.0 


Q 












G30 


11.2 


497. 


48.0 


1.9 


S2 


["} 2MASX source; [*] 2MASSi source; [ c ] SDSS source 





L: LINER, SI, S2: Seyfert 1 & 2, Q: QSO, Int: interaction 
M(H 2 ) and SFE are defined in Sec. 14721 



The CO luminosity for a high-z source is given by 



K km s pc 



where Ico is the intensity in K km s _1 , is the area of the main 
beam in square arcseconds and Di is the luminosity distance in 
Mpc. We assume here that the sources are unresolved, in our 
beam of typically 50-100kpc. We then compute H2 masses using 
Mh 2 = a L'co, with a = 0.8 M (K km s 1 pc 2 ) 1 , for ULIRGs. 
The molecular gas masses are listed in Table[3] Although it could 
be advocated that a different conversion factor should apply to 
some of the galaxies, we always refer to the M(H2) mass, di- 
rectly proportional to CO luminosity, for the sake of compari- 
son. For those galaxies where we observed two CO transitions, 
we used the CO luminosity of the lower transition to calculate H2 
masses. In our sample, most galaxies have CO(1-0) data, except 
three galaxies, G19, G23 and G27, which have been detected in 
CO(2-l). We assume that the brightness temperatures are similar 
in the two lines, as expected for an optically thick, and thermally 
excited medium. The problem is more severe for high-z objects, 
where the CO excitation is not well-known. The average CO lu- 
minosity for the 18 galaxies detected is L'co = 2. 10 10 K km s _1 
pc 2 , corresponding to an average H2 mass of 1.6 10 1() M . 

The star formation efficiency (SFE), also listed in Table [3] is 
defined as Lfir/M(H2) in L /M . Since the SFR is related to the 
FIR luminosity as SFR= L FrR /(5.8 1O 9 L ) (e.g. Kennicutt 1998), 
the gas consumption time-scale can be derived by r — 5.8/SFE 
Gyr. 



4.3. Molecular gas excitation 

Three sources have been detected in both the CO(1-0) and CO(3- 
2) lines, and three have upper limits, as listed in Table [4] For the 
data points, we took the peak flux S v of the lines, since the CO(3- 
2) and CO(1-0) have sometimes different measured linewidths, 
which could be due partly to the noise. The corresponding ra- 
tio between the peak brightness temperatures are also displayed 
in Table [4] to compare more easily with the predictions of the 
model. 

The peak flux ratio between the two lines S32/S10, and equiv- 
alently the peak brightness temperature ratio, is a good indicator 
of the average density of the emitting medium, since density is 
the main factor determining the excitation. Another factor is the 
kinetic temperature, which could be linked to the dust tempera- 
ture (e.g. Weiss etal. 2003). In Table [3] we have computed the 
dust temperature deduced from the far-infrared fluxes, assuming 
k v oc 1/, where k v is the mass opacity of the dust at frequency 
v, and p — 1.5. The average dust temperature for our sample 
is 46 ± 5 K. This is comparable to recent results for starburst 
galaxies, which have dust temperatures » 40 K (e.g. Sanders & 
Mirabel 1996, Elbaz et al. 2010). 

If the gas is predominantly heated by collisions with the dust, 
the gas temperature is expected to be lower than the dust temper- 
ature, at low density (e.g. Spitzer 1978). Alternatively, if the gas 
is heated directly from UV photons near star forming regions, or 
by shocks due to turbulence or perturbed dynamics, then the gas 
could be at much higher kinetic temperatures. Since we observe 
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very low excitation temperatures, we consider it unlikely that, in 
average over the beam, the hot molecular gas is dominating the 
emission. We then assume a gas kinetic temperature lower than 
the dust temperature, in the following modeling. 

Using the Radex code (van der Tak etal. 2007), we have 
computed the predicted main beam temperature ratio between 
the CO(3-2) and CO(l-O) lines, for several kinetic temperatures 
and as a function of H2 densities and CO column densities. 
Figure[5]shows these predictions for T^=45 and 20K. The black 
contours delineate the range of observed values. In Table [4] we 
list the derived values for the n(H2) densities, for two values of 
the kinetic temperatures (45 and 20K), and for a fixed column 
density per velocity width. 

We adopted a column density of N(CO)/AV of 7xl0 16 
cm" 2 /(kms _1 ), which assumes that all CO lines are optically 
thick. This number is at the right order of magnitude given the 
high molecular gas masses derived in Table [2] For M(H2) = 
3 10 10 Mq, a typical CO abundance of CO/H 2 = 10~ 4 , and a 
linewidth of 300 kms , this column density corresponds to a 
homogenous disk of 3 kpc in size. Either the emitting CO gas is 
more concentrated, as in nuclear starbursts, i.e. the CO column 
density would be higher, or the gas extent is larger than 3kpc, 
in which case we would have to take the clumping factor into 
acount. In any case, it is likely that the CO lines are optically 
thick. 

The excitation of the CO gas appears quite low, implying a 
low average H2 density in our galaxies. Comparing to the differ- 
ent excitation patterns observed in other high-z starburst galaxies 
(Weiss et al. 2007), our galaxies are among the lowest excitation, 
comparable to the Milky Way or even lower. 

For the estimations, we have assumed that our sources are 
unresolved. Increasing the size of the molecular disks, here sup- 
posed to be point like (wrt to the beam sizes) to for example >6" 
would increase the T/, ratios by >50%, and this would raise the 
required densities by, at most, a similar factor according to Fig. 

H 

It should be kept in mind that error bars are large on the 
observed ratios. The average ratio is lower than 1.8+0.6, taking 
into account the upper limits. However the conclusion of a rather 
low average H2 density is rather robust with respect to the error 
bars, since the predicted flux ratio is increasing very quickly with 
density in the model. We note that even with the extreme hypoth- 
esis of optically thin gas, the observations are only compatible 
with n(H2) < 3 10 3 cirT 3 , since the predicted ratio increases even 
more with density than in the thick case. It is also possible that 
some of the gas has a kinetic temperature much higher than the 
dust temperature, but then the derived H2 density is even lower. 

4.4. Variation with redshift 

Does the molecular gas content of galaxies evolve with redshift? 
It is interesting to compare the CO luminosity of our sample with 
the wealth of data reported in the literature. Figure [6] shows all 
CO measurements as a function of redshift. This figure reveals 
that indeed our sample (full black circles) is filling in the CO 
redshift desert, although not completely. The rise of the CO lu- 
minosity that is observed at high redshift (z >1) in fact begins as 
soon as z=0.2-0.3. 

This variation is meaningful in the sense that only the bright- 
est objects have been selected here. Most of the variation with z 
comes from the fact that there are no extremely luminous objects 
at z< 0.2. In itself, it is already an interesting evolution, that has 
been discussed in previous works at high redshift (i.e. Solomon 
& vanden Bout 2005, Tacconi et al. 2010). The present work ex- 



tends this variation in the intermediate redshift range, and sug- 
gests that the increase in gas content with z might begin as soon 
as z=0.3. The possibility of undiscovered large CO luminosity 
objects locally is not high, given the good correlation between 
CO and FIR luminosity. These objects should have been discov- 
ered as ULIRGs. 

To interpret this evolution, caveats have to be kept in mind. 
At high redshift, at least for some of the sources, the CO lumi- 
nosity could be over-estimated by poorly known amplification 
factors due to lensing. The luminosities have been corrected for 
amplification, when known, but these uncertainties contribute 
to the large scatter. This is not the case for the sample from 
Daddi etal. (2010, green circles) or the sample from Genzel 
et al. (2010, blue asterisks). Another uncertainty comes from the 
CO gas excitation. The CO luminosities of high redshift ULIRGs 
come from the measured high-J lines, and the low-J lines are of- 
ten not known. They could underestimate the H2 mass, since 
the subthermal excitation is likely to reduce their luminosity 
with respect to the local objects, observed in CO(1-0). We think, 
however, that the steep rise at 0.2 < z < 0.6 of the most lumi- 
nous galaxies discussed in this paper does not suffer from these 
caveats (they are detected in majority in CO(1-0) and are not 
lensed). It is interesting to note that the trend seen in Fig|6]is es- 
sentially dominated by three galaxies, Gl, G15 and G18, which 
are particularly strong in CO-emission. No such extreme L'co 
has been found in the local universe. Two of these galaxies (Gl 
and G18) have been shown in Sect. 14.31 to be subthermally ex- 
cited, and are likely extended starbursts with relatively low effi- 
ciency, as presented by Daddi etal. (2008). In paper I, we also 
derived an extended gas disk for G4 with the Plateau de Bure 
observations. Table [^confirms that Gl and G18 have among the 
lowest SFE of the sample. It is conceivable that the conversion 
factor between L'co and M(H2) could be higher in these sources, 
leading to higher gas mass. 

4.5. Correlation between FIR and CO luminosities 

Figure [7] shows the well studied correlation between FIR and 
CO luminosities (e.g. Young & Scoville 1991). The correla- 
tion is non-linear, with the ultra-luminous objects displaying a 
higher FIR luminosity (i.e. star formation) for the amount of gas 
present, as implied by the lines plotted in the Figure. Our sam- 
ple galaxies fit perfectly in this picture, being all above the curve 
Lfir/M(H2)= 100 L /M (corresponding to a consumption time- 
scale of t =58 Myr). One of the two galaxies detected above 
Lfir/M(H 2 )=1000 L q /M (G19) has indications of nuclear ac- 
tivity (see Type in Table [3]), but the second (G27) has none. We 
note that they are two of the three galaxies detected in CO(2-l), 
and not in CO(1-0); we have derived their H2 masses with the 
assumption of equal CO luminosity between these two first lines. 
If their gas was sub-thermally excited, their H2 mass could then 
be slightly under-estimated. 

Assuming a dust temperature and the observed 100 /jm 
flux S100, we can derive the dust mass as 

M ^ = 4 -8xlO- 11 ( T ^^M 

= 5(1 +z)- (4+ «Sioo^ pc {exp(144(l +z)/T d ) - 1} M 

where S vo is the observed FIR flux measured in Jy, £>^ pc is the 
luminosity distance in Mpc, B vr is the Planck function at the 
rest frequency vr = vo(l + z), and we use a mass opacity coef- 
ficient of 25 cm 2 g -1 at rest-frame 100 fim, (Hildebrand 1983, 
Dunne et al. 2000, Draine 2003), with a frequency dependence 
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Fig. 5. Peak Tb ratio between the CO(3-2) and CO(l-O) lines versus the H2 density, and the CO column density per unit velocity 
width (Nco/AV) for two values of the kinetic temperature: T& = 45K, the dust temperature (left), andT^ = 20K (right). The black 
contours are underlining the values obtained in the data. The predictions come from the LVG hypothesis in the Radex code. 
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Fig. 6. Measured CO luminosities, corrected for amplification 
when known, but not for gas excitation, as a function of redshift. 
We compare our points (full black circles, and arrows as up- 
per limits) with a compilation of high-z molecular gas surveys, 
and local ones: green triangles are from Gao & Solomon (2004), 
blue squares from Solomon etal. (1997), open circles from 
Chung etal. (2009), blue diamonds from Geach etal. (2009), 
black crosses from Iono et al. (2009), red stars, from Greve et al. 
(2005), green full circles from Daddi et al. (2010), blue asterisks 
from Genzel etal. (2010), and blue full circles from Solomon 
& vanden Bout (2005). For illustration purposes only, the red 
curve is the power law for D.uJQ.ui proposed by Obreschkow & 
Rawlings (2009). 



of /3=1.5. Estimated dust masses are displayed in Table[3] If we 
adopt the low conversion factor of a = 0.8 M (K kms~' pc 2 )~\ 




Log LCO (K.km/s.pc 2 ) 

Fig. 7. Correlation between FIR and CO luminosities, for our 
sample (full black circles, and arrows for upper limits) and the 
other points from the literature (same symbols as in Fig|6]l. The 
3 lines are for L FIR /M(H 2 )=10, 100 and 1000 L /M from bot- 
tom to top, assuming a conversion factor a = 0.8 M (K km s~' 
pc 2 ) -1 . The three lines correspond to gas depletion time-scales 
of 580, 58 and 5.8 Myr respectively. 



the average gas-to-dust mass ratio is 96 for all the detected galax- 
ies. The gas to dust mass ratio could increase to up to 550 if the 
standard (MW) conversion factor is used. For local ULIRGs, this 
number is around 100 (Solomon etal. 1997), and 700 for nor- 
mal galaxies (Wiklind etal. 1995), when calculated from IRAS 
fluxes. 
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We should note that the dust temperature has been measured 
from 60 and 100/mi fluxes, which correspond to (1+z) higher 
frequencies in the rest frame of the galaxies. Therefore we are 
not sensitive to the very cold dust (~ 10K). 

Table 4. CO gas excitation 



S32/S10 
[peak] 



Tb32/Tbl0 

[peak] 



n(H 2 ) n(H 2 ) 
[cm 3 ] [cm 4 ] 
45K 20K 



Gl 


2.2 ± 1.1 


0.25 ±0.1 


63. 


200. 


G3 


2.6 + 0.8 


0.29 + 0.09 


79. 


250. 


G17 


1.9 + 0.5 


0.21 ± 0.05 


46. 


140. 


G18 


< 1.3 


<0.14 


<20. 


<63. 


G22 


< 1.0 


< 0.1 1 


<12. 


<40. 



n(H 2 ), for T t =45K and 20K, and N co /AV=7xl0 16 cm^kms" 1 ) 



4.6. Activity of the galaxies 

We have made a census of the different activities occuring in our 
sample galaxies. The last column of Table [3] indicates nuclear 
AGN activity and/or perturbed morphology. These have been de- 
rived from the SDSS images, some of which are shown in Figure 
[8] We note that among the 12 non-detections, there are only 4 ac- 
tive objects, all being Seyfert 1 or quasars, while among the 18 
detected ones, there are 12 active objects, and most of the time 
they are LINERs, Seyfert 2, and show signs of galaxy interac- 
tions and mergers. 



4.7. Star formation efficiency 

In the following we adopt the definition of SFE = Lfir/M(H2), 
with a constant CO-to-H2 conversion factor. The average SFE 
in our sample is 555 L /M , 3 times higher than that of the lo- 
cal ULIRGs (170 Lq/Mq). It should be kept in mind that some 
galaxies could have a different conversion factor, and this uncer- 
tainty affects our conclusions. Also, it is possible that the local 
SFR tracers evolve with time, and that, for a given SFR, differ- 
ent amounts of gas are consumed in star formation, if the IMF 
(Initial Mass Function) is different for earlier and younger galax- 
ies. However, no strong evidence has been found until now for 
a significantly changing IMF, and the star formation laws are 
remarkably constant over redshift, as discussed by Genzel et al. 
(2010) and Daddi etal. (2010). We plot the SFE versus L' co in 
Figure [9] versus Lfir in Figure [TT] and versus the dust tempera- 
ture in Figure [12] 

What is obvious in all these figures is that galaxies of our 
sample are among the most efficient forming stars, and G19 and 
G27 even lie above the starbursts at high redshift. They are not 
among the most gas-rich, according to the CO luminosity. They 
could be experiencing a burst due to galaxy interactions. This 
is the case for G19 (Figure |S). No image is available for G27 
(which lies outside of the footprint of the SDSS). Note that, as 
expected, there is a much larger correlation between SFE and 
Lfir, than with L'co- The presence of large amounts of gas is not 
a sufficient condition to trigger a starburst, and another hidden 
factor is the extent of the spatial distribution of the molecular 
gas. 

To estimate the relative gas fractions of our sample galaxies, 
an estimation of their stellar and/or dynamical mass is required. 
We have tried to estimate the stellar mass from optical and 
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Fig. 9. Star formation efficiency SFE=L F ir/M(H 2 ), versus CO 
luminosity, assuming the same CO-to-H2 conversion factor a = 
0.8 M (K kms -1 pc 2 )~'. All symbols are as defined in Fig [6] 
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Fig. 10. The star formation rate (SFR) obtained from the far in- 
frared luminosity, versus the stellar mass of galaxies in our sam- 
ple (full black circles), compared to the sample of Da Cunha et 
al (2010, full red squares) and Fiolet et al (2009, full green tri- 
angles). 



near-infrared luminosities, taken from the literature, mainly the 
SDSS and 2MASS catalogs. The multi-wavelength luminosities 
were K-corrected according to the colors (e.g. Chilingarian et al. 
2010), and stellar masses estimated according also to the colors 
(Bell et al 2003). Stellar masses were found between 10 1() and 
10 12 M Q or somewhat higher in the case of quasars. Figure [TOl 
displays the star formation rate, derived from the infrared lumi- 
nosity, versus the stellar mass, in comparison to the ULIRG sam- 
ples of Da Cunha et al (2010) and Fiolet et al (2009). Our sample 
points follow the general trend, with the bias of strong SFR, due 
to our selection on L(FIR). The gas fractions derived from these 
stellar masses show large variations, with 3 objects around 1 %, 
but most of them are between 10% and 65%. Another estimation 
is the dynamical mass, derived from the observed CO linewidths. 
This can only be a rough estimate, since neither the inclinations 
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Fig. 8. Optical color images from the Sloan Digital Sky Survey (SDSS, http://www.sdss.org/) of 15 of our sources. They are all 
detections, except G10 and G21. Each panel is 50"x50" in size, and is centered on the galaxy coordinates of Table Q] 



of the galaxies, nor the extent of their CO emission, are known. 
Adopting a typical radius of 3kpc, the dynamical masses are in 
the range of 10 11 M Q , and the derived gas fraction are in gen- 
eral a few percent, with a large scatter, up to 60%. Note that the 
derived stellar masses are on average larger than the dynamical 
masses; this is due to the radius we have selected (3kpc) to es- 
timate the dynamical masses. For the same velocity width, the 
dynamical mass scales as the radius. Firm conclusions cannot 
be drawn until CO maps are available. Indeed, it is possible that 
the sample has a wide range of radial extents and, consequently, 
conversion factors. 

Finally, there is a correlation between the dust temperature 
and the SFE in Figure [T2j it is conceivable that a concentrated 
starburst heats more efficiently the dust around it. However, the 
correlation becomes more scattered at higher redshifts, where 
the largest efficiencies occur. It is possible that our estimation 
of the dust temperature is not as accurate for these more red- 
shifted objects. It is interesting to note the evolution of SFE with 
redshift, where our two most efficient starbursts are clearly no- 
ticeable, as shown in Figure Qj] We qualitatively compare this 
evolution with the cosmic star formation history, as compiled by 
Hopkins & Beacom (2006), from different works in the litera- 
ture, and complemented at very high redshift by the gamma-ray 
burst (GRB) data of Kistler et al. (2009) and the optical data 
(Lyman-Break Galaxies, LBG) from Bouwens et al. (2008). The 
SFE logarithmic variations should be a combination between 
variations of the SFR and of the gas fraction in galaxies. It is 
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Fig. 11. Same as Figure [9] but versus Lfir. 



interesting to superpose the observed logarithmic curve of these 
SFR variations to have an indication of the relative role of the 
various parameters. The schematic curve in log reproduces the 
relative variations, whatever the vertical units, and can be arbi- 
trarily translated vertically. Our points correspond to the most 
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drastic change in this curve, and our following study at 0.6 < z 
< 1 should give more insight in this epoch. 
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Fig. 12. Same as Figure|9] but versus Tj, for the sources where it 
could be defined. For the submillimeter galaxies, the red star cor- 
responds to the averaged SFE, with the mean dust temperature 
of 35.5 found by Kovacs etal. (2010). 
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Fig. 13. Same as Figure [9] but versus redshift. The red curve 
is a schematic line summarizing the cosmic star formation his- 
tory, from the compilation by Hopkins & Beacom (2006), com- 
plemented with the GRB data by Kistler et al. (2009) and the 
optical data from Bouwens et al. (2008). The red curve is loga- 
rithmic and only indicative of relative variations of the star for- 
mation rate per cubic Mpc as a function of redshift, and can be 
translated vertically. 



5. Discussion and conclusions 

We have presented our search for CO-line emission in a sam- 
ple of 30 ULIRGs, selected between 0.2 < z < 0.6 to fill the 



gap or "CO redshift desert" between z=0.1 and 1. We intend 
to cover the second part 0.6 < z < 1.0 in a following work. 
Our detection rate is ~ 60%. We find that some of the galax- 
ies possess large amounts of molecular gas, much larger than 
local ULIRGs. Considering the evolution with cosmic time, it 
appears that the huge amounts of gas, common at high redshift, 
begin to disappear at z~ 0.3. This drop in gas content is reminis- 
cent of the drop in the star formation history of the universe, 
which may imply that the change in star formation is due to 
a change in gas content. There are good reasons to think that 
galaxies are more gas rich at high redshift, and also that their 
gaseous medium is denser. The sizes of galaxies are predicted to 
vary as (l+z)~', and the implied higher gas pressure could in- 
crease the H2/H1 ratio. Following semi-analytical simulations, 
Obreschkow etal. (2009) followed the H2/H1 ratio statistically 
over 30 million galaxies, and its cosmic decline has been mod- 
elled as Qh 2 /^hi °c (1 +z) L6 by Obreschkow & Rawlings (2009). 
This law appears to reproduce grossly the decline in the maxi- 
mum L'co with time, as shown in Figure [6] Some galaxies of 
our sample, however, lie significantly above this envelope. The 
increase with z in the H2 content of galaxies might occur already 
at lower z than this model predicts. 

Five of our galaxies were observed in both CO(3-2) and 
CO(1-0) lines, allowing an estimation of the excitation temper- 
ature. They appear all very low, similar to what is observed in 
the Milky Way, or more normal galaxies, but also some local 
ULIRGs (Radford etal. 1991). These ULIRGs could be similar 
to those discovered by Daddi etal. (2008) at redshift z~1.5. If 
a galactic conversion factor was adopted for these galaxies, as 
suggested by Daddi et al. (2010), their H2 mass would be even 
higher, and they would stand out even more in the cosmic H2 
abundance. 

We emphasize that the choice of the conversion factor is cru- 
cial for the interpretation of the results. We have adopted the 
ULIRG value proposed by Solomon etal. (1997), and derive 
large SFE and short consumption time-scales for the gas. These 
SFE values would be lower if a higher, i.e. the Galactic, con- 
version factor is used. However, for consistency with previous 
studies we prefer to assume only a single value for the conver- 
sion factor. The latter could vary with the extent of the molecular 
component. We do not yet have spatial information on the CO 
emission and future interferometer observations are required to 
constrain the conversion factor further. 

We have compared the star formation history and the red- 
shift evolution of the SFE. It is expected that the latter evolves 
as a combination of the SFR and gas fraction evolution. It is 
likely that the star formation decline between z=l and z=0 is 
partly due to the declining star forming efficiency. For galaxies 
of our sample, the star formation efficiency (SFE) appears very 
high, in comparison to the most active starbursts at different red- 
shifts. This supports a high contribution of the SFE to the star 
formation variations with redshift, although we are observing an 
increase of efficiency of the most extreme objects. It is possible 
to compare the observed time gradients in the cosmic star forma- 
tion rate, and those in the extreme SFE in FigureQj] We observe 
a significant gradient in SFE, but however less steep than in the 
star formation history. The latter requires also a strong variation 
in gas content. 

The very efficient star forming objects (ULIRGs) dominate 
the star formation at high redshift (z < 1 .5, Lefloc'h et al. 2005), 
and less extreme objects (LIRGs) dominate later on (Caputi et al. 
2006), which might explain the strong decline in efficiency be- 
tween z=l and 0. It appears that the range of redshift studied 
here is just where the most massive objects continue to form 
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stars with unprecedented efficiency, before the sudden drop due 
to star formation quenching (e.g. Springel et al. 2005). 
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